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We describe various dynamical processes aimed to create and fuse Majorana fermions in an in- 
homogeneous Kitaev's wire. We show that, with the undesired excitations suppressed by inhomo- 
geneity, fusion of Majorana fermions will result in universal measurable excitations, which can serve 
as a direct verification of Majorana fermions fusion rule. Moreover, we design a protocol to detect 
the oscillatory tunneling between two unpaired Majorana fermions. Most remarkably, our proposal 
is valid for transverse Ising chain, which provides a promising way to test nontrivial properties of 
Majorana fermions experimentally. 
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Majorana fermions (MFs) are self-conjugate quasipar- 
ticles (7T = 7) [1], and non-Abelian anyons obeying ex- 
otic braiding statistics and fusion rule [2,3]. Recent years 
have seen much excitement over MFs, not only because 
of their peculiar properties, but also due to possible ap- 
plications for topological quantum computation [4] . Cre- 
ating, manipulating and detecting MFs experimentally 
remains a great challenge, although many theoretical 
schemes for that have been proposed [5-15]. As far as 
their realization is concerned, an encouraging progress 
has been made for one dimensional systems, especially 
for semiconducting wires [16], where a zero- bias conduc- 
tance peak (ZBCP) [17] and a fractional Josephson ef- 
fect [10, 18] have been recently measured [19, 20]. How- 
ever, it still remains controversial whether the ZBCP 
arises from the MFs [21] and the fractional Josephson 
effect can also be observed in non-topological supercon- 
ductor junctions [22]. Hence up to now we lack such evi- 
dence for MFs, which would unambiguously and straight- 
forwardly demonstrate their fusion rules and braiding 
statistics. 

Experiments in quantum simulation offer a promising 
platform to study and mimic many physical systems from 
condensed matter to cosmology [23] . The theory suggests 
that with a spin-orbital coupling in ultracold atoms, one 
can realize topological superfiuids [8] as well as Kitaev's 
wire [24]. However, these theoretical proposals are still 
hard to implement experimentally. On the other hand, 
transverse Ising chain model (TIM) presents itself as a 
system which is easy to realize and to control in quan- 
tum simulation experiments. It has been realized in a 
large variety of experimental systems, including ultracold 
atoms [25], trapped ions [26], and Josephson circuits [27]. 
TIM is related to Kitaev's wire by a non-local Jordan- 
Wigner transformation. The danger is, however, that 
such non-local transformation may render MFs obscure. 

In this Letter, we show that one can verify MFs fusion 
rule by measuring non-equilibrium dynamics of an inho- 



mogeneous TIM chain. In particular, we discuss non- 
equilibrium dynamics in Kitaev's wire, involving creat- 
ing, manipulating and fusing MFs. We introduce an in- 
homogeneity in the system, which helps to suppress un- 
desired excitations and to manipulate the MFs. We high- 
light that the fusion of MFs would give rise to a promising 
result in TIM chain, an emergence of universal number 
of excitations. The universal emergent excitations may 
serve as a touchstone of nontrivial fusion rule of MFs. 
We also find that oscillatory tunneling between unpaired 
MFs may show up and can be measured via oscillatory 
number of emergent excitations. 

The Hamiltonian and its features.- In the following, 
we consider an open TIM chain with an inhomogeneous 
magnetic field. The Hamiltonian can be expressed as 
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The Jordan- Wigner transformation: of = 2qc; — 1 and 

of = n}=i a j( c l + c i); ma P s TIM into Kitaev's p-wave 
superconducting wire model: 
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H = -5>i(2ctc,-l)-X; J( Cj -ct)(c i+ i+ct +1 ) (2) 



where Cj denotes spinlcss fermion operator, hj can be 
considered as a site dependent chemical potential. 

For TIM, the ground state is doubly-degenerate fer- 
romagnetic state when \hi\ <C J. For the opposite 
limit \hi\ 3> J, the chain will be in the spin polarized 
state with all spins oriented along the magnetic field. 
The two phases are separated by the quantum critical 
point (QCP); for the homogeneous case hj = h, QCP 
is at \h\ = J. Correspondingly, when \h\ < J, the 
Kitaev's wire is in a topological phase, which has two 
zero energy MFs located on the both ends. The exis- 
tence of MFs can be easily understood by considering 



2 



the limiting case h = 0. Expressing the fermion op- 
erator Cj = (72,7-1 + *72j)/2 with MF operator, which 
satisfies 7] = ~f s and {^ s ,lt} = 25 st , one obtains H = 
—iJ~^2j = x l2jl2j+\- Recombining the MF operator into 
fermion operator Cj — (-f 2 j + 172^+1 )/2, the Hamilto- 
nian becomes H = — J ^2^=^ (2cjCj — 1). Therefore, the 
groundstate of the chain is a quantum state with all Cj 
fermions occupied, and with two unpaired MFs 71 and 
72JV localized at both ends. For two unpaired MFs, they 
obey the MFs fusion rule, fusing into a superposition of 
the vacuum state and one fermion c [4]: 



7x7 = + . 



(3) 



Verifying the fusion rule of MFs comprises two steps: 1) 
creating two unpaired MFs, 2) measuring the fused state. 
Even though MFs are created in pairs, one can still have 
two MFs from different pairs, which are unpaired and 
satisfy the MFs fusion rule in Eq. (3). 

Inhomogeneous potential.- For a topological phase 
with well defined MFs, two fused MFs may enter either 
into vacuum or into one fermion state. Since both are 
the groundstates of the system, it is hard to distinguish 
them. To measure the fused state more transparently, 
one may drive the system across the QCP to a non- 
topological phase, where the one fermion state becomes 
an excitation, making it easier to be differentiated from 
the vacuum state. However, the vanishing energy gap 
and divergent relaxation time at QCP will inevitably cre- 
ate many excitations, whose density is determined by the 
Kibble-Zurek mechanism (KZM) [28, 29]. For instance, 
imposing a linear quench h = — J + t/tQ from t < to 
t > 0, will produce excitations whose number n by the 

^ /2 

usual KZ mechanism scales as n a NtQ 1 [29]. These 
excitations will obscure the MFs fused excitation. 

In order to measure the fusion excitation of MFs more 
effectively, one needs to suppress the KZM excitations. 
An intuitive way is to do the experiment in a system of 
small size. For a system with size L, the minimum gap 
is proportional to 1/L, and KZM excitations may vanish 
if the quench process is sufficiently slow. However, the 
smaller the system the greater is the overlap between the 
MFs located at its opposite ends, making their measure- 
ment impossible. An alternative way to suppress KZM 
excitations is to consider an inhomogeneous system [30]. 
This physical result can be understood qualitatively as 
follows: when an inhomogeneous system undergoes a 
quench process, the critical point will be crossed locally 
and the MF associated with boundary between the topo- 
logical and non-topological phases will shift as the quench 
process goes on. It is important that the energy spectrum 
of the bulk of the inhomogeneous system always has fi- 
nite gap during the quench process. Consequently, the 
inhomogeneous setup not only suppresses the bulk KZM 
excitations but also helps one to manipulate the MFs. 

Numerical simulations.- First we consider two simple 



quench processes, where the chemical potential is inho- 
mogeneous and varies with time: 

Process I: h n (t) = a 2 n 2 + h + t/t Q (4) 
Process II: h n {t) = a 2 {n - N/2 - 1/2) 2 + h + t/t Q 

where a denotes the coefficient of parabolic confining po- 
tential. Here we choose the system size N = 100. With 
an increase of a, the minimum gap during the whole adi- 
abatic quench process also increases, making it easier to 
suppress the KZM excitations. The term t/tQ in (4) rep- 
resents the quench, with £q being the quench time, which 
determines the rate of change of the chemical potential. 
During the two quench processes, the chemical potential 
will be ramped from h n < — J to h n > J. Using the 
Bogoliubov transformation = J2n( u n™ c i + v nm,c n ), 
we can diagonalize the Hamiltonian for any given time, 
and find that for both quench processes there is always 
an energy gap [31]. 

To study the non-equilibrium dynamics, we consider 
the BCS groundstate of the chain \ip) = fj a\ |0) /N , and 
the state after the evolution is \ipf) = Yl a \ it /)|0) / 'N 
with 



al(t f ) = Ualu\ U 



T jexp 



H(t)dt 



(5) 



where Nq is normalization constant, U (t) is time evolu- 
tion operator, T is time ordering operator. We also diag- 
onalize the final Hamiltonian H(tf) = '}2, m {E m b\ yl b m — 
Embmbln), with E m < 0. Then the number of excitations 
can be written as 



(0) 



Fig. l(a)-(b) illustrate the number of excitations for 
two quench processes, respectively. It is clear that KZM 
excitations will be greatly suppressed for long quench 
time tQ. In particular, process I gives rise to no exci- 
tations while process II yields excitations with universal 
number 1. To further understand the above numerical 
results, we should carefully consider the two quench pro- 
cesses. In quench process I, the chain is prepared in a 
topologically trivial phase (Fig. l(c)-I). Then by grad- 
ually increasing the chemical potential, we create two 
MFs by pulling them out of the vacuum at the right end 
of the system (Fig. l(c)-II). As the quench goes on, the 
two MFs will separate initially (Fig. l(c)-III), then ap- 
proach each other (Fig. l(c)-IV), and finally fuse at the 
left end (Fig. l(c)-V). The whole process only involves 
two paired MFs, no matter how they move they will al- 
ways fuse into the vacuum, yielding no excitations. This 
process can't be used to detect MFs fusion rule, but it 
can serve as an adiabatic way to evolve the whole chain 
into different quantum phase. 

In the quench process II four MFs participate. As de- 
picted in Fig. l(d)-II, two pairs of MFs, (71,72) and 
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FIG. 1. The number of excitations after the quench process I 
(a) and II (b). A cartoon picture for the two quench processes 
I (c) and II (d). Topological phase regions are marked in 
light cyan. MFs with the same color, (71,72) and (73,74), are 
paired MFs, which will fuse into the vacuum state. Here we 
choose J = 1. 



(73,74), are pulled out of the vacuum on both ends. 
When the potential is further increased, two unpaired 
MFs, 72 and 73, may fuse at the potential center (Fig. 
1 (d)-III to IV), into a superposition of the vacuum |0) 
and one fermion |1), as shown in Eq. (3). Then the 
quantum state can be written as (|0)|0') + |l)|l'))/\/2, 
where |0) and |1) refers to the fusion state of 71 and 74 
to be vacuum or a fermion. When the whole chain is 
in the topological phase (Fig. l(d)-IV), only |1') is an 
excited state, then n ex = 0.5. Similarly, when the whole 
chain is ramped non-topological (Fig. 1(d)- VI), both |1) 
and |1') are excited states, thus n ex = 1. This univer- 
sal single excitation is a direct consequence of the MFs 
fusion rule. 

Measuring MFs fusion rule in TIM. - The groundstate 
of TIM is a BCS type state with all eigenstates of the Bo- 
goliubov quasiparticles filled, thus the excitation number 
created by the quench process is the same as Eq. (6). 
However, due to the quantum fluctuations, the excita- 
tion in TIM is usually a highly entangled superposition 
of many spin configurations, which is difficult to detect 
experimentally by local measurements. 

To clarify the excitations in TIM unambiguously, one 
may tune the system parameters to two limiting cases: 
hi = or J = 0. At hi = 0, the groundstate of TIM is 
just an exact ferromagnetic state with spin-spin correla- 
tions of nearest neighbors: ((TjcrJ +1 ) = 1. An excitation, 
the domain wall, is located on the bond + 1) with 
(o-]a] +1 ) = -1. When J = 0, the groundstate of TIM 
is an exact paramagnetic state with spins aligned paral- 
lel to the magnetic field, and the excitation is a flipped 
spin. In the quench process II, after fusing 72 and 73 
(Fig. l(d)-IV), the corresponding TIM chain is in the 



ferromagnetic phase. Then one tunes slowly the mag- 
netic field to 0, and measures the spin-spin correlation 

(Ei«+i>* one ma y obtain (E^Ii 1 ofof+i) = N-2, 
signifying the emergence of 0.5 excitation. Similarly, af- 
ter the fusion of 71 and 74 (Fig. l(d)-VI), the corre- 
sponding TIM chain is in a paramagnetic phase. Then 
one tunes slowly the spin-spin coupling J to 0, and mea- 
sures the magnetization along the z direction, one may 
have (E*Li °f) — N — 2, signifying the emergence of 1 
excitation. Thus, the universal fusion excitation of MFs 
can be measured in TIM. However, the emergent excita- 
tion here is not well localized in real space [31], which is 
inconvenient for experimental detection. 

To facilitate a detection of the MFs fusion rule, we de- 
sign a quench process III (shown schematically in Fig. 
2(a)), which requires more delicate control on potential 
but the emergent excitations in TIM becomes perfectly 
trapped, as illustrated in Fig. 2(b). Initially we prepare a 
topological trivial phase with the potential monotonically 
increasing along the chain (Fig. 2(a)-I), and then slowly 
increase the potential so that the whole chain enters into 
the topological phase with two paired MFs (71,72) ap- 
pearing at the both ends (Fig. 2(a)-II). This process 
must be adiabatic as was discussed for the quench pro- 
cess I. Then we adiabatically change the potential to be 
harmonically distributed (Fig. 2(a)-III). After that, we 
decrease the potential, creating two paired MFs (73,74) 
at the center (Fig. 2(a)-IV). Finally, 71 and 73 fuse at 
the left end, 72 and 74 may fuse at the right end (Fig. 
2(a)-V to VI). This process III is similar as process II, 
where we fuse two unpaired MFs, nevertheless, our nu- 
merical calculation show that in a remarkable way the 
fused excitations (flipped spins) are perfectly trapped at 
the both ends of the system [31]. 



(a) 



I 




II 























Y2 



IV 


Yi 


y 2 


y, 


u 



V 




h 




y. 


Y4 



VI 



(b) 



Quench III 



iiiiiiii+JiiMi* 



FIG. 2. (a) The quench process III. The topological phase 
regions are marked in light cyan. The MFs with same color, 
(71,72) and (73,74), are paired MFs, which will fuse into the 
vacuum state, (b) A schematic state evolution of the quench 
process. Initial state (groundstate): all spins pointing in the 
— z direction. The final state: an equal-weighted superposi- 
tion of the groundstate and excited state with two end spins 
flipped into the z direction. 



Oscillatory tunneling of unpaired MFs.- When two 
MFs (7a, 7b) become close so that the interaction between 
them is significant, H' = —iAjajf,, this can be approxi- 
mated as a time evolution process with a time evolution 
operator U = exp(-iH't). For two paired MFs, denoted 
as 71 — 172, the time evolution only gives a phase factor, 
U("fi — 272)^ = cxp(i2At)(^i — 172)- For two unpaired 
MFs (72, 73), the situation is quite different. We suppose 
that 72 is paired with 71 and 73 is paired with 74. Then 
the state under the time evolution becomes 

[/(71 — 172)^ = 71 — icos(2At)^f2 — isin(2At)j3, 
U{ia - il3)U f = 74 - i cos(2A£)73 + i sin(2Ai)7 2 . (7) 

Therefore, for two unpaired MFs within a finite range, 
there will be oscillatory tunneling between them. 

In fact, we can design a process IV to measure such 
oscillatory tunneling between two unpaired MFs. The 
process IV is similar to the process II plotted in Fig. 
1(d). However, here we only ramp the chain from Fig. 
l(d)-I to Fig. l(d)-III, and pause at Fig. l(d)-III. The 
MFs 72 and 73 are separated by I lattice spacings, having 
a finite coupling — 1^4(^)7273 between them [32]. Then 
we let the system evolve for certain time duration f e , 
such that 72 and 73 will tunnel to each other, with the 
resulting states 72 = cos[2A(Z)i e ]72 + sin[2A(Z)i e ]73 and 
73 = cos[2A(^)i e ]73— sin[2^4(/)i e ]72. After that, we slowly 
drive the chain back from Fig. 1 (d)-III to Fig. l(d)T, 
where we will fuse (71,72) and (73,74) at the left and 
right ends. The fused excitations can be easily measured 
in TIM, since they will be perfectly localized at the edges 
of the chain, similar to Fig. 3(a). One can obtain that 
the number of excitations created by the whole process: 

n ex = 1 — cos[2A{l)t e + tp ] (8) 

The phase factor ipo comes from the fact that, when two 
unpaired MFs approach and separate with each other, 
the tunneling between them yields a nonuniversal phase 
factor (ipo). This phase factor is difficult to determine 
analytically. However, in experiment one can make the 
process (from Fig. I(d)-I to III as well as from Fig. 1(d)- 
III to I) identical, then phase ipo will be the same for 
the quench processes with different time duration t e . In 
Fig. 3(b), we plot numerical results for the excitation 
number created by such process IV, which clearly shows 
sinusoidal oscillations. The oscillation period of n ex is de- 
termined by r ~ Tt/A(l), with which we can calculate the 
two MFs coupling strength A(l), as shown in Fig. 3(c). 
We find that \nA(l) can be well fitted by a quadratic 
function. 

Braiding MFs through T-junctions in TIM.- It has 
been shown that one can use T-junctions in Kitaev's wire 
to implement braiding operations of MFs [3]. However, 
T-junctions design for TIM will fail to detect the braid- 
ing statistics of MFs, since the model with T-junctions is 



4 




te I 



FIG. 3. (a) A schematic state evolution of the quench process 
IV. Initial state (groundstate): all spins pointing in the —z 
direction. The final state: a superposition of the groundstate 
and excited state with two end spins flipped into the z direc- 
tion. Here /3 = At e + ipo/2. (b) The oscillation of excitation 
number as a function of duration time t e . (c) The coupling 
strength A(l) between two unpaired MFs as a function of dis- 
tance /. We take the following parameters, a — 0.02 and 
J = 1. 



not one-dimensional and so the Jordan- Wigner transfor- 
mation does not map it to a fermion model with a local 
Hamiltonian. 

Experimental realization - Up to now, several experi- 
mental systems may become candidates to test our theo- 
retical proposal. One promising approach is to simulate 
TIM using ultracold atoms on optical lattice, where TIM 
with additional magnetic field in the x direction has al- 
ready been constructed [33]. Moreover, the technique 
for single site measurement in ultracold atoms becomes 
available [34]. On the other hand, TIM has been success- 
fully simulated in the trapped ions system, though with 
only a small number (n ~ 8) of spins [26, 35]. A recent 
work [36], which has hundreds of spins with variable- 
range spin-spin Ising-type interactions simulated, pave 
new way on large-scale quantum simulations for TIM. 
Another way to simulate spin chain models is to use 
Josephson superconducting circuits. A decoherence time 
as long as 0.1ms has been achieved in Josephson circuits 
which is sufficient to test our theoretical predictions [37] . 

Summary-We have studied quench dynamics of Ki- 
taev's wire under the inhomogeneous potential. We pro- 
pose and numerically confirm that, with an inhomoge- 
neous structure of the quench, the fusion rule of MFs 
can be unambiguously detected with a remarkable sup- 
pression of the KZM excitations. The fusion of MFs will 
result in universal excitations, which may serve as a di- 
rect evidence of the MFs fusion rule. We also find that 
oscillatory tunneling between two unpaired MFs will cre- 
ate excitations whose number oscillates with the time of 
evolution. We show that in TIM one can measure both 
universal and oscillatory emergent excitations. Our the- 
oretical proposal provides a promising way to verify the 
nontrivial properties of MFs experimentally in ultracold 
atoms, trapped ions or Josephson junctions. 
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Supplementary material 



ENERGY SPECTRUM 



Using a Bogoliubov transformation = ^2 n (i 



iC n ), Kitaev's wire Hamiltonian 



H 



N 

£ 



N-l 



^•(2 C t Cj - 1) - Jfa - 4)( c j+i + 4+i) 



(9) 



can be diagonalized into H = 'Yl im {E m al n a m — E^amO^), where E m < 0. Due to the particle-hole symmetry, the 
spectrum has conjugate pairs with negative/positive energy, and we call quasiparticle, a m quasihole. One can 
obtain the energy spectrum during the two quench processes: 

Process I: h n (t) = a 2 n 2 + ho + t/tq 



(10) 



Process II: h n (t) = a 2 {n- N/2 - 1/2) 2 + h + t/t Q 



The energy spectrum is shown in Fig. 4. One can observe that the there is always an energy gap during the quench 
process, and Majorana zero modes will appear when part of the chain is engineered into topological phase. Especially, 
for the process II, four Majorana fermions (MFs) will appear. 




FIG. 4. Energy spectrum of quench process I (a) and II (b). In (a), we take the parameters as a = 2 7 , ho = —1.8, N = 100, 



and t/tQ from to 3. In (b), we take the parameters as a = 2 
we only plot several states with energy around 0. 



-1.7, N = 100, and t/tQ form to 3. For simplicity, 



NUMBER OF EMERGENT EXCITATIONS 



The eigenstate of the system can be constructed by N particles, each of them is chosen from (aj„,a m ) conjugate 
pairs, which reads, 



H') = u^ n °iio>/jvo 

jGA i<£N~A 



(11) 



where Nq is the normalization constant. To calculate the number of excitations created by the quench process, 
we start with a groundstate, \ipo) = Ili a !|0)/-^0: evolve the state into |V>/)> with time-evolution operator U = 



\1>f) = W = U~[[al\0)/No = H(Ua\U^\0) 



(12) 



Therefore, we can simply calculate the time-evolution of each quasiparticle a] n: o^itf) = Ua\ n U^ . Apparently, the 
evolved quasiparticle aj„(i/) still satisfy the fermionic commutation relation {a m (tf), a n (tf)} = {(^(tf), a„(t/)} = 0, 
and {aln (tf),a n (tf)} = S mn . 



7 



To implement the time-evolution numerically, one may discretize the time-evolution operator by using the time- 
slicing procedure: 

U = T jexp -i jf ' H(t)dt j « cxp [-iff(*i)Ai] (13) 

with At -C 1. It should be noted that it is crucial to retain the unitarity of exp [— i?J(tj)At] throughout the calculation: 

cxp [-iif^At] = Aexp(-iAAt)A* (14) 

where H{U) = AAA*, A is a unitary matrix AA^ = I and A is a diagonal matrix. 

Meanwhile, one can diagonalize the final Hamiltonian H(tf), obtains Bogoliubov quasiparticles (&{,••• ,b N ), and 
Bogoliubov quasiholes (b%, ■ ■ ■ ,bisr). The groundstate of the final Hamiltonian is composed with all the quasiparticles 
(b\, • • • , fojy). Therefore, the excitation number in the evolved state \ipf) is n ex = 'Ylnty f\bib\\ip f) . We may rewrite b n 
as superposition of aj n (t/) and a m (tf): 

(t f )+V nm a m (tf)]. (15) 

m 

The excitation number can be written as 

Ei 7 ?™! 2 ( 16 ) 



n.m 



EXCITATIONS IN TIM 



To measure the fused excitations in TIM unambiguously, we may either tune the magnetic field hi or the spin-spin 



coupling J to 0. One can measure the spin-spin correlation (afa. 



i+ll 



\){cj+i + cj +1 )) at hi = to study 



the domain wall excitations, as shown in Fig. 5(a). Similarly, one can measure the magnetization (erf) = (2cjcj — 1) 
when J = to witness the flipped spin excitations, as shown in Fig. 5(b)-(c). It clearly shows for the quench process 
II, either the domain wall excitations (Fig. 5(a)) or the flipped spin excitations (Fig. 5(b)) are not localized well. On 
the contrary, for the quench process III introduced in the paper, where the MFs are fused at the edge of the chain 
and may result in the prefect trapping of the fused excitation, as plotted in Fig. 5(c). 
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FIG. 5. (a)-(b) Spatial distribution of the domain wall and flipped spin excitations in quench process II. (c) Spatial distribution 
of the flipped spin excitations in quench process III. 



TUNNELING OF UNPAIRED MFS 



Wc numerically calculate the cases where two unpaired MFs are approaching, and study how the wave packets 
behave. More specifically, as proposed in this paper, we use quench process II to create two pairs of MFs (71,72) 
and (73,74), and bring two unpaired MFs (72,73) close to each other at the center. Then we pause the changing of 
potential, and let it evolve under that situation for certain time duration t e . The numerical results are shown in Fig. 
6. It clearly shows that finite tunneling appears between the two unpaired MFs (72,73) and the tunneling exhibits 
oscillating behaviors: 



72 = cos[2A(Z)t e ]72 + sin[2/l(Z)ie]73 , 73 = cos[2^(/)i e ]73 - sin[2 A(l)t e ]j 2 



(17) 
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where A(l) is the coupling strength of 72 and 73. 

It can also be observed that the MFs (71,74) at both ends are very localized, and shows no time evolution. All 
these facts suggest that the interaction terms like —17172 or —17374 are negligibly small which we can safely ignore. 
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FIG. 6. Time evolution of the wave packets of four MFs. Here we choose a — 0.02, J — 1, N = 100 and 72 is separated 
by / = 35 lattice spacing with 73. In (a)-(d), we choose the time duration t e to be 100, 300, 671, 1586. Here the period of 
oscillation is tt/A(1) w 1486. 



